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Aims. To describe the evaporation status of the extrasolar planets, we propose to consider an energy diagram in which the 
potential energy of the planets is plotted versus the energy received by the upper atmosphere. 

Methods. Here we present a basic method to estimate these quantities. For the potential energy, we include the modification of 
the gravity field by the tidal forces from the parent stars. 

Results. This description allows a quick estimate of both the escape rate of the atmospheric gas and the lifetime of a planet 
against the evaporation process. In the energy diagram, we find an evaporation-forbidden region in which a gaseous planet 
would evaporate in less than 5 billion years. With their observed characteristics, all extrasolar planets are found outside this 
evaporation-forbidden region. The escape rates are estimated to be in the range 10^ gs"^ to 10^^gs"\ with few cases above 
10^^ g s~^ . The estimated escape rate for HD 209458 b is found to be consistent with the lower limit of 10^" g obtained from 
interpretation of the Hi Lyman-a observations. 

Finally, this diagram suggests possibilities for the nature of the recently discovered Neptune-mass planets. We find that 
GJ436b, 55 Cue e and HD 69830b cannot be low mass gaseous planets. With density necessarily above 0.5 gem"'' to survive 
evaporation, these planets must contain a large fraction of solid/liquid material. Concerning GJ876d, we find that it must 
have a density larger than ~3gcm~^ to survive the strong EUV energy fiux from its nearby parent star. GJ 876 d must contain 
a large fraction of massive elements. 



Key words. Stars: planetary systems 



1. Introduction 



Among the nearby two hundreds extrasolar planets iden- 
tified today, there is a particular class of planets very 
close to their parent stars with semi major axis SSO.l AU. 
Among them, there are massive hot-Jupiters (also named 
"Pegasides", Guillot et al. 1996) and small mass planets 
(see for instance, Santos et al. 2004; Rivera et al. 2005). 
These planets constitute about 20% and 5% respectively, 
of the planets identified today. Since the discovery of the 
first of them, 51Pegb (Mayor & Queloz 1995), the issue 
of their evaporation status has been raised (Burrows & 
Lunine 1995; Guillot et al. 1996). Since these planets ex- 
ist and orbit stars which are not particularly young, it 
is clear that the evaporation rate of the observed mas- 
sive hot-Jupiters should be modest enough to not dra- 
matically impact their evolution. At least they can sur- 
vive billion years at the place where they are observed. It 
is extremely unlikely that we are witnessing a transient 
phase in the evolution of extrasolar planets (Baraffe et al. 
2004) . In brief, the discoveries of a large number of massive 



hot-Jupiters^ and very-hot- Jupiters^ led to the conclusion 
that the evaporation of massive planets has to be modest. 

In that frame the discovery that the transiting extra- 
solar planet HD 209458 b is indeed losing mass was mostly 
unexpected (Vidal-Madjar et al. 2003). From the observa- 
tion of 15±4% absorption in Lyman-a during the transit, 
a lower limit of the H i escape rate was estimated to be of 
the order of ~ lO^^gs-^ (Vidal-Madjar et al. 2003; Vidal- 
Madjar & Lecavelier des Etangs 2004). This mass loss rate 
corresponds to a loss of ~0.1% of the total mass of the 
planet over 10 billion years (10 Gyr). Models of the time 
variation of this escape rate lead to slightly larger values 
for the total fraction of mass loss over the whole stellar 
life (see e.g., Lecavelier des Etangs et al. 2004; Baraffe et 
al. 2004). 

Different detailed models have been developed to un- 
derstand the observed evaporation process and to estimate 
the escape rate (e.g., Yelle 2004; Tian et al. 2005, Garcia 
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^ usually defined as massive planets orbiting in less than 
10 days 

^ usually defined as planets orbiting in less than 3 days {e.g., 
Gaudi et al. 2005) 
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Muiioz 2006), to be compared with the observational con- 
strains (Vidal-Madjar et ah 2003, 2004). Finally, it has 
been suggested that evaporation can lead to significant 
modification of the planets nature in the case of planets 
with a mass only a fraction of the mass of the observed hot- 
Jupiter and orbiting close to their parent stars (Lecavelier 
des Etangs et al. 2004; Baraffe et al. 2005). 

However, models were applied to specific cases; it is not 
yet possible to conclude on the order of magnitude of the 
possible evaporation of planets which are being discovered 
without the effort to model in detail each particular case. 
The purpose of the present paper is not to review previous 
works and extended modeling efforts. Here we propose an 
alternative approach to obtain general characteristics from 
basic observed parameters. This approach allow statistical 
conclusion without the need to model in detail every in- 
dividual planet for which detailed physical characteristics 
are not available. We propose to outline the main char- 
acteristics of a given planet that allow concluding on the 
general trend of its evaporation status, escape rate and 
lifetime related to the erosion through evaporation. 

Here we present the idea of an energy diagram to ad- 
dress these issues (Sect. |21), and the way to simply esti- 
mate the energetic characteristics of a given planet (Sect.O 
and^J. Then an evaluation of the lifetime and escape rates 
is presented in Sect.[S| The resulting diagram is shown in 
Sect. A diagram corresponding to the escape rate is 
given in Sect.[71 Uncertainties and evolutionary tracks are 
presented in Sect. |S| and El Finally we conclude on the 
possible nature of the recently discovered Neptune mass 
planets (Sect. [TU|) . 

2. Potential energy of planets versus stellar 
energy flux 

A plot of the lifetime of planets as function of their masses 
and orbital distances has been calculated by Lecavelier 
des Etangs et al. (2004). The purpose of this plot is to 
show that the identified planets are stable against evap- 
oration. For the existing planets, the lifetime, defined as 
the time to escape the full mass of the planet, is always 
longer than ^ 10^" years (Fig. 5 of Lecavelier des Etangs 
et al. 2004). Attempts have been made to extrapolate from 
this plot possible escape rate in the case of planets discov- 
ered later (HD 149026b, Valenti 2005). It would be ex- 
tremely valuable to have a quick diagnostics on the evap- 
oration rate of a given planet to predict if there are any 
chances to observe and constrain the evaporation through, 
e.g., Lyman-a observation during a transit (Lecavelier des 
Etangs et al., in preparation). However, in the lifetime 
plot of Lecavelier des Etangs et al. (2004), the energy flux 
from the star is assumed to be the one of HD 209458 b. 
The energy flux from other stars is certainly different de- 
pending on the stellar type, rotation period, etc. The sec- 
ond critical parameter in this plot is the assumed plane- 
tary radii. In the lifetime plot of Lecavelier des Etangs et 
al. (2004), the radius of HD 209458b (1.35 Jupiter radii) 
was used to normalize the scaling laws for the planetary 



radii because at the time of this work HD 209458 b was 
the only hot-Jupiter whose radius had been measured. 
However, the radius of HD 209258 b is now considered as 
anomalously large (Guillot 2005; Lecavelier des Etangs 
& Vidal-Madjar 2006). From observations of ten transit- 
ing planets (Bouchy et al. 2005b, McCuUogh et al. 2006) 
and detailed modeling of the planetary structure (Guillot 
2005), 1.1 Jupiter radii is certainly more appropriate for a 
1 Jupiter mass planet close to its parent star. This makes 
HD 209458b with its 1.35 Jupiter radii and 0.69 Jupiter 
mass clearly a peculiar object. In the following we will 
therefore use a fit to the Guillot's model. A new plot has 
been calculated with this new information (Lecavelier des 
Etangs & Vidal-Madjar 2006). However, the planetary ra- 
dius depends also on parameters other than the planetary 
mass and orbital distance, for instance, the stellar type 
and the corresponding effective temperature of the cen- 
tral star and its age. It is thus impossible to uniquely 
estimate lifetime and escape rate from the position of a 
given planet in the orbital distance versus mass diagram. 
As a consequence, this plot has to be revisited. 

In summary, since each planet has different radius and 
different energy flux entering onto the upper atmosphere 
from stars of different type, the orbital distance versus 
mass diagram is not appropriate to conclude on the es- 
cape rate and life time of a variety of planets. We have 
to define the parameters (or combination of parameters) 
which allow us to conclude on that matter. 

After the realization that HD 209458 b is losing mass, 
many modeling efforts have been pursued to better under- 
stand the escape in the particular case of hot-Jupiters. The 
common results of these numerous models is that almost 
100% of the EUV/Lyman-a energy flux which is converted 
into heat is then transferred into escape (see, for instance, 
Yelle 2004 and Tian et al. 2005). Even magnetic fields 
which can inhibit ion escape are found to have modest 
impact. Yelle (2004) found that the net escape rate with 
ion escape inhibited by magnetic field is only 30% smaller 
than in its reference model. Similarly, effects of dissocia- 
tion and ionization of H2 are found to be negligible. The 
resulting escape rates of these various models are different 
not because these models assumed different mechanisms 
but because of the different illumination, planetary cross 
section for the energy flux, and fraction of this energy 
converted into heat to the upper atmosphere. 

As a result, the evaporation status of a set of different 
planets orbiting different stars can be obtained by a com- 
parison of the energy deposited in the upper atmosphere 
and used to escape the planet gravity, i.e. to compensate 
for the (negative) potential energy of the atmospheric gas, 
with this potential energy. We therefore propose to calcu- 
late the position of the extrasolar planets in a diagram 
displaying the potential energy versus the stellar energy 
deposition. We now evaluate these two quantities in the 
two next sections to plot the corresponding diagrams in 
Sect, m and □ 
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3. The extreme ultraviolet illumination from the 
central stars 

The present and past illumination from the parent stars 
is a key parameter which needs to be accounted for. First, 
as the planet-hosting stars are cooler than early F stars in 
the main sequence, the stellar energy emitted in the near 
UV range is negligible (e.g., Kurucz 1993). Moreover, near 
UV and optical energy is absorbed deep in the bottom and 
dense atmosphere. On the contrary, neutral components 
of atmospheres have large cross section in the extreme 
ultraviolet (100A< A <1200A). Photons in the extreme 
ultraviolet wavelength range (EUV) are thus easily ab- 
sorbed by the tenuous upper atmospheres of planets. For 
instance, in the terrestrial planets of the Solar system, the 
EUV energy flux is the main driver of the escape (e.g.. 
Chamberlain & Huntcn 1987). Therefore, for simplicity 
and in agreement with the result of all detailed models of 
the evaporation of hot-Jupiters, we can assume that the 
energy flux deposited in the upper atmosphere and trans- 
ferred into escape energy mainly comes from the EUV. 

However, because of the interstellar absorption, there 
are very few observational constraints on the EUV flux 
from the stars in the solar neighborhood. On the other 
hand, although stellar EUV fluxes are related to stellar 
rotational velocities (Wood et al. 1994), most of the stars 
harboring extrasolar planets have no published rotational 
velocities or only upper limits. In most cases, the infor- 
mation on the planet-hosting stars is limited to the stellar 
type. We have hence chosen to consider the very extreme 
ultraviolet luminosities Ls2 as a function of the stellar 
type, as measured by Rosat in the S2 channel between 
110 and 200 A. We then scale these luminosities to the 
solar EUV flux of ~4.6ergcm~^ s^^ at 1 AU as given by 
Ribas et al. (2005). This solar value is obtained consider- 
ing mid-solar cycle data of 1993 as representative of the 
Sun at the average activity. The EUV flux, i^EUv(l AU), 
received by a planet by unit area at 1 AU is hence calcu- 
lated using 



Feuv(IAU) =4.6 



L 



S2 



erg cm ^ s ^ . 



(1) 



The scale between L52 luminosities and EUV is justified 
by the fact that both the ~ 100 — 200 A and the extreme- 
far UV fluxes are emitted by the same region in the chro- 
mosphere and the transition region where the temperature 
is observed to increase from the photospheric to the coro- 
nal values. 

Hodgkin and Pye (1994) give median S2 luminosities 
as a function of the spectral types, with value of 10^^-^, 

respectively. Concerning the F type stars, most of the F 
stars known to harbor planets are late type F8v or F9v 
stars. Because early F stars remain chromospherically ac- 
tive and fast-rotating during their main sequence phase, 
the accuracy in the measurement of their radial velocity 
is strongly hmited. It is yet impossible to detect planets 
around them (F. Pont, private communication). But these 



early F stars contribute to the mean value of the observed 
S2 luminosities of the whole class of F stars. We therefore 
decided to consider a typical S2 luminosity of 10^^ '^ erg s~^ 
for F6 and F7 stars harboring planets, and an S2 luminos- 
ity of lO^^'^ergs"^ for later FSv and F9v stars which are 
more similar to GOV type stars. 

Using these S2 luminosities and the solar EUV lu- 
minosity, we find Feuv(1 AU) = 14.7, 4.6, 4.6, 14.7 and 
2.9ergcm-2s-i for F6-F7, F8-F9, G, K, and M stars, re- 
spectively. Except if explicitly mentioned (Sect|Hl), these 
values will be used throughout the paper. These fiuxes can 
be converted to the energy per unit of time, dE-Euv / dt, 
received by a planet of radius Rp with an orbital distance 
Up using: 

dEE^y/dt = ttRI (^) Feuv(1 AU). (2) 



4. Potential energy 



To evaluate the total lifetime of a planet in an energetic 
point of view, we need to consider the total potential en- 
ergy, Ep, of the planet considered as a whole, from its sur- 
face to its deep interior. If this energy is given is en?, the 
ratio with the received energy expressed in erg per billion 
of years will provide an estimate of the total time needed 
to cancel the (negative) potential energy, in other words, 
the evaporation lifetime. On the other side, to evaluate 
the escape of the atmospheric gas at the present time, we 
need to consider the potential energy of the gas only in 
the atmosphere and per unit of mass. The ratio of this 
quantity with the energy received will provide the mass 
escape rate. We see that we have two different quantities 
which need to be estimated: the potential energy of the 
whole planet and the potential energy per unit of mass in 
the atmosphere 

4.1. The potential energy of the whole planet 

The potential energy of a gaseous sphere of radius Rp is 
given by 

m(r) dm{r) 



G 



dr 



-dr 



(3) 



where G is the gravitational constant and m(r) is the 
mass inside a radius r. For a uniform density and a mass 
Mp, we obtain the classical formula Ep = —3GMp/5Rp. 
But in general, the potential energy depends on the inter- 
nal structure. For a planet with a mass Mp and a radius 
Rp, an analytical approximation can be found assuming a 
polytropic approximation of the equation of state for the 
planet interior. Using a polytropic index n = 1 (P cx p^, 
de Pater & Lissauer 2001) the potential energy is found 
to be: 



Ep = --GMl/Rp 



(4) 

The details of the calculation are given in Appendix 
Numerically this converts into: 

2 / n \ -1 



E„ 



^2.9 X 10 



43 ( \ ( 



R3 



up 



erg. 



(5) 
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4.2. The potential energy per unit of mass in tine 
atmosphere 

Another interesting quantity is the potential energy of the 
atmospheric gas per unit of mass. This potential energy 
per unit of mass is simply given by (see Appendix : 



dEp(^i^tm) / dm = -GMp/Rp. 
Numerically this converts into: 



dE. 



'p(atm) 



/dm 



1.9 X 10' 



M 



Jup 



Rn 



ergg 



(6) 



(7) 



We note that the potential energy per unit of mass 
of the atmospheric gas is much higher than its kinetic or 
thermal energy, even if the gas is escaping at high veloc- 
ity. For instance, even for gas escaping at high velocity of 
lOkms"^, the kinetic energy is 5 x 10^' ergg~^, which has 
to be compared with the typical value of —dEp(^^^_^-^/ dm ~ 
2^Qi3gj.gg-i ^Yie known extrasolar planets (Sect. d). 
Even for very-hot-Jupiters in which tidal forces indeed 
play a significant role fSect. W^ . the potential energy is 
expected to be always larger than the kinetic energy. The 
kinetic energy of the escaping gas is thus neglected in the 
present analysis. 

4.3. Potential energy and tidal forces 

In the case of hot-Jupiters and very-hot-Jupiters, the grav- 
ity from the nearby parent star modifies the geometry 
of the potential energy. In those cases, the iso-potentials 
around the planet at high altitude are no more spherical 
and are open to the infinity when the potential reaches 
the value of the Roche lobe. Therefore the energy needed 
to escape the planet can be defined by the energy needed 
to reach this Roche lobe. In the case of Mp <C M*, the 
potential energy of the atmospheric gas per unit of mass 
is accordingly modified by 

dE'p(^t,^) /dm = dEp(^^^^) /dm + Atidaid-Ep(atm) /dm, (8) 

where (see Appendix Ib)| 

34/3 gmI^^'m^/'' 



atm) 

/dr. 



(9) 



Op being the orbital distance of the planet and Af* the 
mass of the central star. 

The potential energy of the whole planet as defined in 
Sect. 14.11 is also modified by tidal forces. We can define 
Atidai-E'p by 

E'p^Ep + Atidai^p. (10) 

Fortunately, because Eq. El does not depend on the plan- 
etary radius, the change introduced by the tidal forces 
does not depend on the planet internal structure. Indeed, 
we have 



itidal-C/p 



E„ — 



3^/3 GMl^^m,{r)y'^ dm{r) 



dr 



dr. (11) 



It is easy to find that for the total potential energy of a 
planet we have 



Atidal-Ep 



3^/3 GA/y Afp'/' 
10 a„ 



(12) 



Because Atidaic'£'p(atm)/c^'7i is positive, the energy deposi- 
tion needed to escape the planet gravity is decreased by 
the tidal forces. For the identified hot-Jupiters, we find 
that the depth of their potential well is decreased by the 
tidal forces by typically 30 to 50%. For instance, we find 
that Atidaid-Bp(atm)/|di^p(atm)| ~50% for OGLE-TR-56b 
and --33% for HD 209458 b. Keeping the hypothesis that 
the escape energy is mainly used to escape the potential 
well of the planet, this leads to a similar increase of the 
escape rate (Sect El- 



5. Lifetime and escape rate 

Assuming that the energy deposited by the EUV flux is 
used to extract the gas which escapes the potential well 
of the planet with an efficiency e = 1, we can derive an 
escape rate, dm/dt, and a corresponding typical lifetime 
for a given planet. Assuming ^ c'^'p(atm) ^ ^dE-E\]Y, and 
neglecting the tidal forces (d-Ep(atm) ~ dEp^^^^^)), we find: 



dm/dt — 



dE-Evw/dt 



dEp(^^tm) / dm' 
From Eq. El we derive the escape rate 

dTu/dt = Feuv(1 AU) 




(13) 



(14) 



Numerically, this gives 

dm/dt =3.7 xlO'f ^EUv(lAU) 
\ 4.6 erg cm ^s ^ 



Jup 



( Mp 



Jup 



(15) 



1 AU 



Concerning the lifetime, we can similarly assume 
—E'p = e Jp*^ di^EUV to obtain a typical duration, ti, which 
corresponds to the time needed to fill the whole potential 
well of the planet assuming the current energy fiux. This 
duration can hence be called a "lifetime" ; the nature of 
a planet can change through evaporation if ti is shorter 
than the typical lifetime of the parent star. If the energy 
received by the planet, dE-Eijv/dt, is assumed to be con- 
stant over time, from the above equation we derive the 
lifetime ti: 



-K 



dEEvv/dt' 



(16) 



However, when the planetary mass evolves, the radius 
and the corresponding energy flux captured by the planet 
changes, consequently also the real lifetime changes. A 
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more detailed analysis is possible by considering the vari- 
ations of the escape rate as a function of time. This leads 
to a new "lifetime", t2, which can be defined by 

-K = 

(tMj) ' r ""^pitTiFEVvil AV)it'))dt', ^^^^ 

where to is the initial time of the system. At this stage a 
useful quantity is the mean value of the EUV energy flux 
from the central star < dE^w/dt > defined by: 

< dEEuv/dt > (t) = 

/ a, y2 ^ jl7rR,{t'nFEUvilAm'))dt' (18) 

which includes the variations of the EUV flux and the 
variation of the planetary radius which both are larger 
during the early stage of the planetary system (Guinan 
k Ribas 2002; Ribas et al. 2005; Burrows et al. 2000). 
We can define a correction factor for time variation of the 
energy flux received by a planet, 

^ " dEEvv/dt • ^ ' 

Hence we have 

t2-l-Hi. (20) 

We now have to evaluate the impact of this 7 correction 
factor. First, the variation of the EUV flux as a function of 
time has been thouroughly analyzed for G stars by Ribas 
et al. (2005) within the "Sun in time" program. A recent 
work reports similar results for K type stars (Lakatos et al. 
2005). Roughly speaking, the EUV flux can be assumed to 
follow a power law in the form F cx t^", where a is close 
to 1. For G type stars, the 1-1200A EUV flux is found to 
have a ~ 1.23, while the Lyman-a has a — 0.72 (Ribas 
et al. 2005). For K type stars, the stellar Civ emission 
lines, which are believed to emerge from the same coronal 
region as the EUV flux, follows a power law with a ^ 0.94. 
We therefore can consider a power law with a = 1 as a 
good approximation (see also Vardavas 2005). Assuming 
that the increase of the UV luminosities in the past has 
to be limited and should "saturate" , we consider that the 
luminosities do not vary at ages below 100 million years 
(Lakatos et al. 2005). With these values, and neglecting 
for the moment the variation of the planetary radius, we 
find that for a = 1.23 we have 7 = 6.5 at i = 5 Gyr, and 
for a = 0.75 we have 7 = 4. We can therefore consider a 
power law with a = 1 as a good approximation. 

The variation of the planetary radius is a more complex 
matter because this depends on the planetary character- 
istics like mass and orbital distance. Nonetheless, empir- 
ically we found that the published time variation of the 
planetary radius (Burrows et al. 2000, Guillot 2005) can 
be well fitted by a simple scaling law in the form : 

i?p(t) = i?oo+/3i"°'. (21) 



We find good fits with /3 = 0.2R^ Gyr°-^ for Mp Z 
0.3 Mjup, and with /3 = 0.3i?oo Gyr°-^ for Mp « 0.1 Mj^p. 

Using a = 1 for the power low of the luminosity varia- 
tions and P — 0.2i?oo Gyr" "^ for the variation of the plan- 
etary radius, we find 7 = 6.8 at t = 2 Gyr (the typical age 
of the F6-F7 stars harboring planets, Saffe et al. 2005) and 
7 = 8.0 at t = 5 Gyr (taken as the typical age of later type 
stars). Using for less massive planets (3 = 0.3i?oo Gyr°'^, 
we find 7 = 8.7 at t = 2 Gyr and 7 = 10.0 at t = 5 Gyr. In 
the following we consequently decided to consider 7 = 7 
for F6-F7 stars, and 7 = 8 for later type stars. 

If we consider that the EUV luminosities can be ex- 
trapolated with the function in the past earlier than 
at the age of 100 million years, down to the age of 30 
million years, 7 is then increased by about 25%. We thus 
conclude that the correction factor for time variation of 
the energy flux does not strongly depend on the epoch of 
this "saturation" of the EUV luminosities in the past. 

6. The energy diagram 

With the equations derived in the previous sections, we 
can now plot diagrams of the potential energy versus en- 
ergy flux (Fig.^. In particular Eas.l^ and 1121 are used to 
derived the potential energy, E^, and Eq. ^| (or a com- 
bination of Eqs. 121 and I19|l is used to derive the mean 
EUV energy received per unit of time by the planets, 
< dEEvv/dt >. We used the characteristics of the 182 
identified planets as tabulated by Schneider (2006) on 
June 15*'\ 2006. 

The only remaining difficulty in the drawing of the 
plot lies in the mass determination. Except for the case of 
transiting planets (Sect. IS)), most of the planets have only a 
determination of Mp sini, where i is the orbit inclination. 
We compensate for this effect by multiplying the Mp sin i 
by a factor \/2. This \/2 factor is taken because the real 
mass Mp is statistically smaller than ^/2Mp sin i in 50% 
of the cases, and larger in 50% of the cases. The result is 
thus a spread of the position of the planets in the energy 
diagram. For transiting planets because sin i « 1, we keep 
the measured mass without applying a correction factor. 

From these planetary masses and using the tabulated 
types of the parent stars we estimate the planetary ra- 
dius by fitting the radius versus mass curves calculated by 
Guillot (2005). This model assumes a central rock core and 
a gaseous envelope following an hydrogen-helium equation 
of state with an helium mixing ratio Y = 0.30. The esti- 
mated radii should not be sensitive to the evaporation, 
because they result from a calculation of the planet inte- 
rior at several bars pressure not affected by the evapora- 
tion operating much higher in the atmosphere; this is con- 
firmed by the absence of correlation of the observed radii 
of transiting planets with age and EUV flux received dur- 
ing the planets' lifetime (Guillot et al. 2006). An indeter- 
niination remains for the very low mass planets (Neptune 
mass planets) because their radius (or density) as a func- 
tion of their mass strongly depends on their nature. As a 
first approach and to plot the energy diagram in Fig. 
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Fig. 1. Plot of the potential energy of the extrasolar plan- 
ets as a function of the mean EUV energy received per 
billion of years, < dE-Euv / dt >. To keep the planets with 
the smallest orbital distances in the left part of the dia- 
gram, the direction of the abscissa axis is chosen with the 
largest value of the mean energy flux toward the left. The 
182 identified planets are plotted with symbols depending 
on the type of the central star: triangles for F stars, filled 
circles for G stars, diamonds for K stars and squares for 
M stars; planets orbiting class III stars are plotted with 
empty circles. The absence of planets in the bottom right 
region can be explained by the detection limit of the ra- 
dial velocity surveys. There are indeed very few planets 
with radial velocities of the parent stars below the lines 
for which the radial velocity is 10ms~^. The radial ve- 
locity iso-curves at 10 ms^^ have been plotted for various 
stellar types (see Appendix [nj ; the iso-curves for F8-F9, 
G and K-type stars are found to be nearly superimposed 
and are plotted with a single line. 

The striking result is the absence of planets in the bottom 
left region which corresponds to light planets (small —E'p) 
at short orbital distances (large < dEEw/dt >). A plot 
of the lifetime line at i2=5Gyr, shows that there are no 
planets in this part of the diagram simply because this 
is an evaporation-forbidden region. Planets in this region 
would receive more EUV energy than needed to fill the 
potential well of the planet, and evaporate in less than 
5Gyr. 



we consider that these planets have a typical density of 
2 g cm~'^ for 15 Earth mass and 6 g cm"-^ for 6 Earth mass. 
A detailed analysis of these planets in this diagram will 
be presented in Sect. EH 

The first plot (Fig^l shows the total potential en- 
ergy of the planet as a function of the mean energy flux 
< dE^uv /dt > . We decided to plot the mean energy flux 
on the abscissa axis in reverse order with larger values in 
the left. This choice is done to keep the planets with short 
orbital distances in the left part of the diagram for consis- 
tency with previous plots like, for instance, the diagram 
of Lecavelier des Etangs et al. (2004). The opposite of the 



potential energy, —E'p, is larger for more massive planets, 
and massive planets are therefore found in the upper part 
of the diagram. In some way this diagram is close to the 
classical diagram of mass versus orbital distance. However, 
in the present diagram the position reflects the impact of 
the evaporation on the planetary nature. 

We plot the diagram with the 182 identified extrasolar 
planets. The striking result is the absence of planets in two 
different regions of the diagram. The absence of planets in 
the bottom right region can be explained by the detection 
limit of the radial velocity surveys. We plotted the limit 
for radial velocity of 10ms~^. As expected, we see that 
very few planets have radial velocities of the parent stars 
below this limit. 

The most important region devoid of planets in the 
diagram is the bottom left region. This region should con- 
tain light planets at short orbital distances. However a 
plot of the lifetime line at t2=5Gyr, shows that there are 
no planets in this part of the diagram simply because this 
is an evaporation-forbidden region. Planets in this region 
would receive more EUV energy than needed to fill the 
potential well of the planet and evaporate the gas, leaving 
a remaining core, an evaporation remnant (also named a 
"chthonian" planet; Lecavelier des Etangs et al. 2004). 

From the position in the diagram, the typical lifetime 
of a given planet can be rapidly extracted. If the mean 
energy fiux KdE^uy / dt> is given in unit of erg per billion 
years, and the potential energy is given in unit of erg, the 
simple ratio of both quantities provides the correspond- 
ing lifetime in billion of years. In the diagram, lifetime 
isochrones are straight lines. In the Fig. ^ we plotted the 
line corresponding to the lifetime of 5 Gyr. 

7. Diagram of the potential energy of the 

atmosphere and the corresponding escape rate 

With the above equations (Eqs.[21Eland|51) we can also plot 
the potential energy per unit mass of the atmosphere as a 
function of the present time energy flux (Fig. [SJ . We again 
superimposed the position of the 182 identified planets. 

Although apparently similar to the previous diagram, 
this new diagram allows the estimate of the actual evap- 
oration rate for a given planet, assuming that most of 
the energy in the extreme ultraviolet is used to escape 
the combined effect of planet gravity and tidal forces. In 
Fig. ^ we considered the total potential energy providing 
a hint on the total escape mass and time. Here in Fig. 12 we 
consider the potential energy per unit of mass (ergg~^) of 
the atmospheric gas versus the present time input energy 
(ergs~^). The ratio of both quantities provides a hint of 
the present escape rate (gs^^). 

We see in the plot that the extrasolar planets must 
have escape rate in the range 10^ to lO^^gs^^. The 
largest escape rates are obtained for 0GLE-TR-113b, 
HD 189733b, HD46375b, HD63454b, and TrES-1, for 
which we obtain escape rates of 3.5 x 10^^, 3.2 x 10^^, 
2.8 X 10", 2.5 X 10", and 2.2 x 10"gs-\ respectively. 
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Fig. 2. Plot of the potential energy per unit of mass of 
the atmosphere of the extrasolar planets as a function of 
the EUV input flux. Unlike Fig. ^ the potential energy is 
per unit of mass at one planetary radius in ergg~^; and 
the EUV input flux is the value at the present time in 
ergs^^. This allows to directly estimate the correspond- 
ing escape rate at the present time in gs~^. The escape 
rates are drawn with straight lines for 10^, 10*, 10^° and 
10^^ gs^^. The 182 identified planets are plotted with the 
same symbols as in Fig. ^ 

These five stars are cataloged as K type stars and there- 
fore assumed to have large extreme ultraviolet luminosi- 
ties. The next largest escape rate is found for an M type 
star, GJ 436 b, for which we estimate a maximum escape 
rate of the order of 2.2 x 10^^ gs~^. 

For HD 209458 b, we obtain an escape rate of 9 x 
]^QiOgg-i^ which is consistent with the lower limit of ^ 
j^Qio gg-i obtained from the interpretation of the Lyman- 
a transit light curve (Vidal-Madjar et al. 2003; Vidal- 
Madjar & Lecavelier des Etangs 2004). 

8. Uncertainties 

For the observed planets, the position in the previous di- 
agrams are obtained by simple estimates of the EUV lu- 
minosities and planetary radii. These rough estimates are 
subject to uncertainties. These uncertainties add some dis- 
persion in the clouds of points of a given diagram; however 
with a limited dispersion, this does not severely affect the 
main trends and the corresponding statistical conclusions, 
like the empty region observed for lifetime below 5 Gyr. 
More critically, when we consider a particular planet, we 
need to have an estimate of the uncertainties introduced 
by above approximations. 

To evaluate the impact of these approximations on the 
planetary radius and stellar EUV luminosities, we used the 
transiting planets. For these planets, we have a measured 
radius and also a better estimate of the EUV luminosi- 
ties through the rotational velocity of the central stars 
V sin i (when measured, in few cases only upper limits are 




25.5 25.0 24.5 24.0 23.5 23.0 22.5 
Logio dEEuv/dt (erg s"') 

Fig. 3. Plot of the potential energy per unit of mass of 
the atmosphere of the ten transiting planets as a function 
of the EUV input flux. For each planet the position is 
calculated using the radius as modeled (diamonds) and 
the radius as effectively measured (circles). For OGLE- 
TR-lOb we used the two possible planetary radii quoted by 
Santos et al. 2006. When available, the EUV luminosities 
estimated from the stellar rotational velocity are plotted 
with triangles. 

available). Because the rotational axis of the stars is ex- 
pected to be close to the orbital axis of the planets (see 
for instance Queloz et al. 2000), and for transiting plan- 
ets sini ^ 1, the measured wsini arc in fact true rota- 
tional velocities from which better estimates of the EUV 
luminosities can be obtained. We plotted these transiting 
planets in the diagram of the potential energy per unit of 
mass versus EUV energy flux (Fig. OJ. 

Concerning the radius estimate, OGLE-TR-lOb is a 
special case because this planets has presently two differ- 
ent radius estimates which are significantly different. From 
the analysis of Bouchy et al. 2005a, Santos et al. 2006 ob- 
tain a radius of 1.43 Jupiter radii, while they obtain a 
radius of 1.14 Jupiter radii using the new observational 
data of Holman et al. 2005. Because we cannot discrimi- 
nate between both results, we decided to keep both radii 
in the present discussion. 

We first evaluate the uncertainty on the radius esti- 
mate and its impact on the evaporation rate. Because 
dEpi^gXm) / dm (X Rp^ and dE-E\]\/dt cx R^, the resulting 
escape rate dm/dt is quite sensitive to the planetary ra- 
dius and is proportional to R^. For each planet in Fig. 
we consider the radius obtained from our fit to Guillot's 
models (Guillot 2005; plotted with diamonds) and the ra- 
dius as effectively measured (plotted with filled circles). In 
the ten transiting planets the difference between the two 
corresponding escape rates has a mean value of 0.13 dex 
and a dispersion of 0.22 dex. The largest differences are 
obtained for OGLE-TR-10 (increase of the escape rate by 
0.45 dex with the radius given by Bouchy et al. 2005a), 
HD 209458 b (increase of the escape rate by 0.33 dex with 
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the real radius) and HD 149026 b (decrease of the escape 
rate by 0.41 dex with the real radius). This is due to 
the unexplained excess of planetary radius of OGLE-TR- 
IO and HD 209458 b (Guillot et al. 2005, Laughlin et al. 
2005) and the small radius of HD 189733 b explained by 
its large metallicity and its suspected massive core (Sato 
et al. 2005). To conclude, although the escape rate is quite 
sensitive to the radius estimate calculated from the planet 
characteristics, the transiting planets allows us to estimate 
that this leads to a typical uncertainty of less than 0.3 dex 
in the calculated escape rate. 

The impact of the EUV energy flux estimate is more 
difficult to constrain as we do not have real EUV flux mea- 
sured for any of the stars harboring planets. However, we 
can use another estimate of the EUV luminosities through 
true rotational velocity when available in the published lit- 
erature {i.e. for six of the ten transiting planets). For these 
six planets we calculated the EUV luminosities from their 
rotational velocity using the results of Wood et al. (1994). 
Using Rosat S2 measurements of stars within 10 par- 
sees Wood et al. (1994) found a close correlation between 
the S2 luminosity, Ls2, and the stellar rotational veloc- 
ity, Vy-ot- Ls2 oc v^^^. From this equation, and assuming 
Few = 4.6ergcm~^ s~^ at lAU and i^rot — 2.0kms~^ 
for the Sun, we find an EUV energy flux of: 

Feuv(IAU) = 4.6 ^£2i_^ ergcm-^s-i (22) 

\2.0kms / 

at a distance of 1 AU from a star with a rotational velocity 

For the six transiting planets with published rotational 
velocity we plot the result in Fig. with triangles. In 
these six planets the difference between the two corre- 
sponding escape rates has a mean value of 0.18 dex and 
a dispersion of 0.60 dex. The biggest differences are ob- 
tained for OGLE-TR-10 (increase of the escape rate by 
0.82 dex with the EUV Imninosity estimated from the ro- 
tational velocity) and TrES-1 (decrease of the escape rate 
by 0.87 dex). This is due to the large rotational velocity 
of OGLE-TR-10 with vaini= 7.7kms-i (Bouchy et al. 
2005a), and the small rotational velocity of TReS-1 with 
vs'mi— 1.08 km s^^ (Laughlin et al. 2005). As a result, we 
see that in two of the six cases, we have a large uncer- 
tainty of the EUV luminosity with consequently a rela- 
tively large uncertainty in the escape rate (~ 0.8 dex), for 
the four other cases the dispersion is reduced to 0.4 dex. In 
general, taken into account the two uncertainties (radius 
and energy flux estimates), we find that the differences 
in the estimated escape rates have a typical dispersion of 
0.7 dex. 

9. Evolution tracks 

The diagram described above aims at describing the pos- 
sibility of the planets evolution on billion years time scale 
through the evaporation process. The position in the di- 
agram can be calculated or evaluated at a given time, in 



particular in the present epoch (assuming an age for the 
parent star). It is thus interesting to see how a planet can 
evolve in such a diagram. An interesting aspect of the dia- 
gram is that the position in the diagram in itself contains 
the information of the evolution timescale. 

We thus calculated evolution tracks for few cases of 
planets for which the initial conditions imply rapid evolu- 
tion through evaporation. This concerns light planets at 
short distances (Saturn mass planets at very-hot- Jupiter 
distances). In calculating these evolution tracks, the crit- 
ical parameter is found to be the planetary radius and 
mass at a given age, both evolving with time. Again we 
assumed that the EUV flux varies as and used our fit 
to the Guillot et al. models and the corresponding time 
evolution described by Eq. [21 The result is plotted in 
Fig.H 
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Fig. 4. Plot of the potential energy of the extrasolar plan- 
ets as a function of the EUV input fiux with evolution- 
ary tracks superimposed. This is an enlargement of Fig^ 
Time evolution of planets at constant orbital distance 
ftp — 0.021 AU and ap — 0.04 AU, corresponding to the 
orbital distances of G J 876 d and two Neptune-mass plan- 
ets (GJ581b and 55Cnce), respectively. GJ581b and 
55 Cnc e are on the track of a planet with an initial mass 
Mp(io)=0.22Mjup. GJ876d is on the track of a planet 
with an initial mass of Mp(io) '~-"0.3 Afjup, showing that it 
can be the remnant of a former hot Saturn-mass planet. 



We calculated the time evolution of planets at orbital 
distance Up — 0.021 AU and ap = 0.04 AU, corresponding 
to the orbital distances of GJ 876 d and two Neptune-mass 
planets (GJ 581 b and 55 Cnc e), respectively. We find that 
GJ 581 b and 55 Cnc e can well be the remnant of a planet 
with an initial mass of 0.22 Mj^p. The evolution track of a 
planet with an orbital distance of 0.021 AU and an initial 
mass of 0.3Mjup shows that GJ876d can be the remnant 
of a planet whose initial mass could be up to the mass of 
Saturn. 
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Fig. 5. Mass distribution of planets discovered through 
radial velocity. 



10. Neptune mass planets in the diagram 

A plot of the mass distribution of the extrasolar planets 
shows that Neptune mass planets play a particular role 
(Fig. |SJ|. Ten planets have been found with mass below 
0.067 Mjup (21 Earth-mass), while no planets have been 
identified with mass between 0.067 Afjup and O.llAfjup 
(35 Earth-mass). The presence of this gap is not related 
to a bias in the radial velocity searches, since more mas- 
sive planets are easier to detect (at the same orbital dis- 
tances). This reveals the different nature of these Neptune 
mass planets orbiting at short orbital distances. But their 
nature is still a matter of debate (Baraffe et al. 2005). In 
particular the question arises if they can be the remnants 
of evaporated more massive planets ( "chthonian planets" ) 
as foreseen in Lecavelier des Etangs et al. (2004) ? Other 
possibilities include gaseous Neptune-like planets, super- 
Earth (Santos et al. 2004) or ocean-planets (Kutchner 
2003, Leger et al. 2004). 

The mass of these objects being measured, the free- 
dom for their position in the energy diagram mainly relies 
on the determination of their radius, or, in other words, of 
their mean density, i.e., their nature. As a first approach, 
in Fig. n we considered that these planets have a typi- 
cal density of 2gcm~'^ for 15 Earth mass and 6gcm~'^ 
for 6 Earth mass. Now we can test the density hypoth- 
esis and try to conclude on the possible nature of these 
planets. We plotted the position of these Neptune mass 
planets in the energy diagram with different hypothesis 
on their density (Fig. EJ. We used mean planetary den- 
sity of Pp=6gcm~'^ for a typical density of refractory-rich 
planets which should describe the chthonian and super- 
Earth planets. A lower density on the order of Pp=2 g cm"'^ 
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Fig. 6. Plot of the potential energy of the Neptune mass 
planet as a function of the EUV flux for various planets' 
density. For GJ581b, GJ436b, HD 69830 b, 55Cnce and 
GJ 876 d, and assuming sin i = V2~^, lifetime shorter than 
5 Gyr are obtained for densities below 0.28, 0.55, 0.56, 
0.69 and 3.1gcm~^, respectively. If sini=l (dotted line), 
the critical (minimum) densities are increased to 0.38, 
0.74, 0.78, 0.93 and 4.2gcm-3, for GJ581b, GJ436b, 
HD 69830 b, 55Cnce, and GJ876 respectively. 



can be considered as more plausible for volatile-rich plan- 
ets describing the ocean planets. For gas-rich planets 
we assumed much lower density of pp=0.2gcm~^ and 
Pp=0.4gcm^^, describing planets which should look more 
like irradiated Neptune-like planets (Baraffe et al. 2006 
find densities of 0.1 to 0.25 gcm^"^ for the final stage of an 
evaporated gaseous planets with a 6 Earth mass core and 
final mass of about one Neptune mass). 

In the energy diagram, HD 190360 c, HD 4308 b, 
HD 160691 d, HD 69830 c and HD 69830 d are all located 
in the possible domain with lifetime above 5 Gyr, inde- 
pendently of their assumed density. Even with very low 
densities down to 0.2 gcm^-^, these planets do not receive 
enough energy from their central star to evaporate signifi- 
cantly. For GJ581b only the extreme case pp=0.2gcm~^ 
is in the evaporation-forbidden region. Owing to their or- 
bital distances between 0.09 and 0.63 AU, it appears that 
the evaporation is not efficient enough to modify the na- 
ture if these planets. These planet are not likely to be the 
evaporation remnants of previously more massive planets 

For G J 436 b, 55Gnce and HD 69830 b, density below 
0.4gcm~^ would place these planets in the evaporation- 
forbidden region of the energy diagram. If these planets 
were formerly low density (gaseous giant) planets, they 
must have lost a large fraction of their upper atmospheric 
gas. We therefore conclude that, with low orbital distances 
in the range 0.03 to 0.08 AU, these planets cannot belong 
to the category of low mass gaseous planet (irradiated 
Neptune-like planets). The position of these planets in 
the energy diagram shows that they must contain a large 
fraction of solid/liquid material in the form of volatile or 
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refractory elements. These planets can be the remnants 
of previously more massive planets (chthonian planets) or 
alternatively volatile-rich planets (ocean planets) or rocky 
planets (super-Earths) . 

It is noteworthy that the conclusion for HD 69830 b is 
consistent with the one reached by Lovis et al. (2006) who 
concluded from formation and migration models that it 
must have a mainly rocky composition. 

For GJ876d, even the density of an ocean planet cre- 
ates some problem. GJ 876 d needs to be a planet dense 
enough to be located above the ^2 — 5Gyr lifetime limit. 
This planet requires a density larger than 3.1gcm~'^ for 
its atmosphere to survive. GJ 876 d could be a big rocky 
planet, like a super-Earth, or a refractory remnant of 
a previous more massive planet (an evaporated ocean 
planet?). However, we note that its parent star is an M 
star, and we have very little information about the EUV 
flux of non active M dwarfs; even if every M star has pe- 
riod of strong activity remains unclear. A more detailed 
analysis of the energy flux from this star and its possible 
consequences on this intriguing planet is needed. 

In brief, the energy diagram allows us to trace 
three different categories for the ten presently identified 
Neptune mass planets. For six of them, the EUV input en- 
ergy seems not strong enough to affect significantly these 
planets; we cannot conclude on their nature. For three 
other planets, it appears that they cannot be a kind of 
low mass gaseous planets. These planets must contain a 
large fraction of sohd/hquid material. Finally, GJ876d 
must be dense enough to survive the strong EUV energy 
flux from its nearby parent star. This planet must contain 
a large fraction of massive elements. 

11. Conclusion 

To describe the evaporation status of the extrasolar plan- 
ets, we proposed to consider an energy diagram in which 
the potential energy of the planets is plotted versus the 
energy received by the upper atmosphere. We presented 
a basic method to estimate these quantities. For the en- 
ergy flux, we proposed a first approach using the type of 
the parent stars or, if available, their rotation period. For 
the potential energy, we included the modification of the 
gravity field by the tidal forces from the parent stars. This 
description allows a quick estimate of both the escape rate 
of the atmospheric gas and the lifetime of a planet against 
the evaporation process. We found that the region of life- 
time below 5 Gyr is devoid of planets, which reinforces the 
belief that this diagram is an appropriate description of 
the evaporation. In the evaporation-forbidden region, the 
evaporation can modify the nature and the characteristics 
of the planets such that they move to region of the di- 
agram corresponding to longer lifetime and consequently 
become long-lived. 

The escape rates estimated from the diagram are found 
to be in the range lO'^gs^^ to lO^^gs"^, with few cases 
above lO^^gs"^. Although this is a crude estimate, the 
estimated escape rate for HD 209458 b is found to be con- 



sistent with the lower limit of 10^" gs^^ obtained from 
interpretation of the Hi Lyman-a observations. 

Finally, this diagram allows the description of the pos- 
sible nature of the recently discovered Neptune-mass plan- 
ets. We found that GJ 436 b, 55 Cnc e and HD 69830 b can- 
not be low mass gaseous planets. These planets must con- 
tain a large fraction of solid/liquid material. Concerning 
GJ876d, we argued that it must be dense enough to sur- 
vive the strong EUV energy fiux from its nearby parent 
star. GJ 876 d must contain a large fraction of massive el- 
ements. 
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Appendix A: The potential energy of a gaseous 
planet 

The potential energy of a gaseous sphere of radius Rp is 
given by 

^0 ' 



dr 



(A.l) 



where G is the gravitational constant and m(r) is the mass 
inside a radius r. The internal structure, characterized by 
m(r), is given by 



dP_ 
dr 



-pir 



Gm{r) 



where p{r) is the density given by 

— — — 47rr p(r), 
dr 



(A.2) 



(A.3) 



and the relation between the pressure P{r) and the den- 
sity p{r) is given by the equation of state. An analytical 
solution can be found if we assume that the equation of 
state follows a polytropic with an index n — 1, which is 
found to be a good approximation for the gaseous plan- 
ets interior (de Pater & Lissauer 2001). In that case, the 
equation of state is 



P = Kp^ 



(A.4) 



where if is a constant. The equations IA.2I IXTSl and lA.41 

can be analytically solved. This gives the density profile 
within the planet (de Pater & Lissauer 2001): 
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the Roche lobe. Along the star-planet axis, the Roche lobe 
is located at a distance r = rpjochcj where rpocho follows: 



En 



(A.7) 



Also using Eqs. EH EH and EH the planet mass, Mp is 
given by 



dr 



(f'Rochc) = 0. 



(B.6) 



Mp = m{Rp) = m{TT/C'K) 



Again assuming that r <^ ap, and using Ea. lB.3l andlB 
we find 



(A., 



/i = 3 



Finally, putting Eq. IA.6I and IA.8I into I A. 71 we find the 
potential energy of a gaseous sphere in which the equation 
of state follows a polytropic with an index n = \: 
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from which we derive 
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Appendix B: Modification of the potential energy 
by tidal forces 

In the field of a planet of mass Mp, at a distance r from 
the planet center, the potential energy field, x, follows the 
equation 



Finally, the new potential energy of a unit mass of the 
atmospheric gas of a planet, dE'^^^^^-^/ dm, defined by the 
opposite of the energy needed to escape the planet gravity, 
is now the energy to reach the Roche lobe: 



X{r) 



GM„ 



(B.l) 



dE'p(^tra)ldm =x'{Rp) - x'lj'Rocho) 

= [x{Rp) - x(?'Rochc)] 

+ [Ax(i?p) - Ax(rRochc)] 



(B.9) 



where G is the gravitational constant and C is a constant. 
The potential energy of a unit mass of the atmospheric gas 
of a planet, dEp^g_im)ldm, can be defined by the opposite 
of the energy needed to escape the planet gravity, in other 
words to reach infinite distance from the planet: 

dEp(^,tm)/dm = x{Rp) - X{^) = -GMp/Rp. (B.2) 



In the case /z <C 1 and Rp <C Op, using Eq. IB. II and IB. 71 
we derive for the first term of Eq. IB. 91 
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where Rp is the radius of the planet. This demonstrates Yov the second term of Eg. iB 
the well known Eq.El given in Sect. 14.21 

However, in a star-planet system, the potential energy 
field is modified by the star gravity and centrifugal forces. 
Along the star-planet axis, the potential energy field, x'j 
follows the equation: 



using Eg. IB. 51 we derive: 
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x'{r) =x(^) + ^X{r) 
= - G{Mp + M^ 
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where r is the distance from the planet center in the direc- 
tion of the star on the star-planet axis, Op is the star-planet 
distance, is the stellar mass and fi is defined by 



Here we implicitly assumed {r^ocho / Rp)^ S> 1, which is 
valid as the smallest values in extreme cases of very-hot- 
Jupiters are always larger than 10. 

As a result, combining Eg. lR!2llR9llBT0l and lBTTl we 
find 



dEp^^tm)/dm 



M„ 



Mp + M^ 

In the case r ap and /i <C 1, we can derive 
Ax(r) 



GiMp + M^) 



(B.4) 

+ c. 

(B.5) 



dEp(^i,tm)/dTn 
9G{Mp + M,) / 

^R,ochc 



2ar, 



(B.12) 



If we define AtidaidEpf^^^^) / dm by 



dEpi^^tm)/dm I 



dE, 



p(atm 



)ldm 



itidal 



dE, 



'p{atni)/dm, 

(B.13) 



using Eg. lB.8l and lB.12l we find 



In the modified gravity field, the energy needed to escape 
the planet can be defined by the energy needed to reach 



AtidalC?-E'p(atm) 



/dm 



34/3G(Mp + A/,)/i2/3 



2ar) 



(B.14) 
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That is 



tidal '^-E'p(atm) 



34/3 gmI'^'mI'^ 

/dm w . 



(B.15) 



This resuh is obtained using the potential description in 
the direction of the star on the star-planet axis, but it is 
valid in all direction because the potential is the same on 
the whole Roche lobe surface which is defined as an iso- 
potential. Therefore the energy needed to reach the Roche 
lobe and to escape the planet is the same in all directions. 

Appendix C: The radial velocity iso-curves in the 
energy diagram 

Most of the known planets have been detected with the 
radial velocity method which is not sensitive to low mass 
planets with large orbital distances. Therefore, if the de- 
tected planets are plotted in the energy diagram (Fig. , 
the absence of planets with small absolute value of the 
potential energy (low mass) receiving a low EUV energy 
(large orbital distance) is due to this bias of the radial 
velocity method. The amplitude of the radial velocity of 
a planet-hosting star, I^ad, is given by 



Kid = 



(C.l) 



where G is the gravitational constant, ap the planet orbital 
distance, Mp and Af* are the mass of the planet and the 
star, respectively. Using Eq. ^ and the fact that 



*,EUV 



where P*,euv is the total power emitted in EUV by the 
parent star, we derive 



En 



rad 1 



'/rf^EuvMy 

V -P*,EUV / 



(C.3) 



This equation does not depend on the planetary radius, 
mass and orbital distance. Therefore, for a given stellar 
type, corresponding to a given stellar mass and total power 
emitted in the EUV, the radial velocity iso-curves must 
follow a straight line in the log-log energy diagram of Ep 
as a function of dEEw /dt. 

Most of the instruments have a detection threshold at 
about 10ms~^. Instruments with larger sensitivity (like 
the Harps spectrograph) are recent and allowed the dis- 
covery of low mass planets only with short orbital periods. 
This explains why the detected planets below the limit at 
10 ms""'^ are in the left part of the energy diagram plotted 
in Fig.n 
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